ftpft  ±1*  8^-7 

ESTIMATION  AND  TESTS-OF-FIT  FOR  THE 
THREE  PARAMETER  WEIBULL  DISTRIBUTION 


Richard  A.  Lockhart 
Michael  A.  Stephens 


TECEN1CAL  REPORT  No.  476 
OCTOBER  27,  1998 


Prepared  Under  Contract 
N00014-92-J-1264  (NR-042-267) 

FOR  THE  OFFICE  OF  NAVAL  RESEARCH 


Professor  Herbert  Solomon,  Project  Director 


Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government. 


Approved  for  public  release;  distribution  unlimited 


DEPARTMENT  OF  STATISTICS 
STANFORD  UNIVERSITY 
STANFORD,  CALIFORNIA  94305-4065 


Estimation  techniques  are  given  for  the  throe*  paramo  tar  Weibull 
distribution,  with  tha  location  (or  origin)  paramatar  unknown,  and  possibly 
also  the  shape  and/or  scale  parameters  unknown.  Tests  of  fit  are  described, 
and  tables  are  given  for  the  SDF  statistics  A2,  W2  and  U2,  to  make  the  teats. 
Several  examples  are  discussed. 


Key  words.  EDF  tests:  empirical  distribution  function: 


goodness - of - fit : 


reliability:  survival  analysis. 

DTK  QUALITY  INSPECTED  8 


Aocsstoo  Far 

NTW  CRAW 
OTIC  TAB 
Unannounced 


Justification 


By _ 

Distribution/ 

Availability  Codes 


Dtst 


Avail  and/or 
Special 


SSTXMATXOB  AMD  TESTS  OF  TIT  FOB.  THK  THREE- FAIAMKI1 
MEISULL  DISTSItUTIOM 


1.  IMTEODUCTIOM 

In  this  article  nr«  given  estimation  procedures  and  tests  of  fit  (based 
on  the  enpirical  distribution  function,  or  EDF),  for  the  three -parasM ter 


Veibull  distribution: 


F(x;a,0,m)  -  1  -  exp[-((x-a)/0)  ], 


x  >  a 


where  0  and  a  are  positive  constants.  When  o  is  known,  the 
distribution  is  called  the  two -parana ter  Veibull  distribution,  and  estimation 


procedures  and  goodness  ef  /it. tests,  are  then  very  straightforward:  the  tests 

;  V/T  r7:.i^iabA  I 

are  referenced  in  'Sectioik4>fepl<^lv~p*re  we  concentrate  on  tests  for  use  when 

’  ,f.;  SAi  OlTd  j 

a  is  not  known,  but  •  •ftft-rnrTi*qr£rr  fros  the  sample,  together  with  ■  and 

L...  I 

0  if  necessbry.-  In  these  «lrme>rsy»s ,  problems  can  sosMtlmes  arise  in 

* 

estinating  tfreparanetsrs.sofn  faction  2,  procedures  are  given  to  obtain 

),:V*  • 

.....  .  _  ^  t 

efficient  estinates  which  can  then  be  used  with  the  goodness -of- fit  tests. 

t 

For  reasons  ‘of  space  many  ddtailaof  both  estimation  and  goodness-of-fit 
procedures  have  been  omitted; j some  j>f  these  are  in  Lockhart  and  Stephens 

(1989).  I -  .  ! - ! _ ! 

It  is  worthwhile  to  observe  that  the  three -parameter  Veibull 


distribution  is  a 


of  a  wider  class,  the  generalized  Extreme-Value  (or 


Jenkins  on)  distribution.  This  distribution  is 


F  (x;a,b,c)  -  1  -  exp 


x  >  a 


The  parameter  b  must  be  positive  while  a  and  c  may  be  any  real  numbers. 
The  three-parameter  Veibull  distribution  is  the  subfamily  of  F*  with 
c  >  0.  The  special  case  c-0  is  the  usual  extreme-value  distribution 


F  (x;a,b,0)  -  1  -  exp[-  exp((x-a)/b) )  , 


•  <  x  <  •  ; 
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ie  ariMt  as  the  lisle  of  the  three -pares*  ter  Uoibull  f sally  (1)  as  a  -»  •. 


Kstlmatlon  rtf  par— tort.  Tho  goodness -of- fie  procedures  depend  on  flree 
esclaecing  ehe  peraaecers  in  (1)  by  an  efflclene  aechod,  for  example,  by 
aaxlaua  likelihood.  If  one  had  an  lnflnicely  large  Veibull  saaple,  such 
esclaaees  could  be  found  straightforwardly.  However,  for  finiee  eaaplee, 
usually  snail,  there  will  be  sons  daea  sees  which  give  no  local  maximum  for 
ehe  likelihood.  If  one  were  willing  Co  fie  the  larger  family  (2),  Milan 
likelihood  estimates  could  be  found,  but  with  a  negative  c;  thus  ehe 
resulting  fitted  distribution  will  not  be  in  the  Veibull  f sally.  However, 
practical  workers  in  aany  fields  do  not  wish  to  broaden  the  class  of 
distributions  beyond  the  Veibull,  to  include  (2);  then,  for  practical 
purposes,  the  best  Veibull  fit  will  be  the  distribution  with  c-0  (the 
nearest  non-negative  c),  that  is  to  say,  the  extreae -value  distribution  (3). 
Note  that  c-0  corresponds  to  a—  in  (1)  (and  a  -»  -•  and  $■*•). 

In  Section  2  we  discuss  estimation,  and  bow  to  recognise  this 
difficult  case.  A  formula  is  given  to  calculate,  from  a  given  data  set,  a 
quantity  A,  which,  when  negative,  indicates  that  no  local  maximal  for  the 
likelihood  can  be  found,  and  distribution  (3)  should  be  fitted  to  the  data. 

Cheng  and  lies  (1990)  have  discussed  estimation  problems  for  three- 
parameter  distributions  embedded  in  a  larger  model  and  in  particular,  problems 
for  the  three-parameter  Veibull  distribution,  and  have  pointed  out  the  connec¬ 
tions  between  (1)  and  (2).  They  offer  tests  for  c*0,  with  the  assumption  that, 
if  this  can  be  accepted,  the  extreme-value  distribution  (3)  will  be  fitted  to 
a  data  set;  if  a  zero  c  is  rejected,  estimation  continues  for  the  Veibull 
distribution.  The  implication  is  that,  where  acceptable,  the  distribution  (3) 
will  be  fitted  in  preference  to  (1).  Ve  are  concerned  with  providing  tests  of 
fit  for  the  distribution  (1),  but  estimation  of  parameters  must  come  first;  we 
concentrate  therefore  on  fitting  the  Veibull  distribution  whenever  this  is 


u 


possible,  and  only  Call  back  on  (3)  abort  wo  aust  assume  ■  -  There  are 
naturally  nany  points  of  contact  with  Chang  and  Ilea  (1990),  as  we  shall 
see  below,  although  these  authors  do  not  discuss  tests  of  fit. 

Apart  from  the  problea  of  deciding  whether  a  is  infinite,  it  is 
well-known  that  there  are  other  probleas  of  estimation  when  a  is  not 
known;  for  example,  the  likelihood  can  be  aade  infinite  if  a  also  is 
unknown,  or  if  known  a  <  1.  An  extensive  literature  exists  on  tikis  type  of 
problea  (for  exsaple  Smith,  1985;  Cheng  and  lies  1987).  Smith  (1985) 
discusses  asymptotic  procedures  in  detail,  and  gives  extensive  theory,  to 
which  we  refer  below;  Smith  and  Naylor  (1987)  coapare  Bayesian  and  maximum 
likelihood  estimators  in  a  case  study.  Of  the  results  which  are  known  and 
proved,  most  of  them  concern  asymptotics,  where  the  sample  sixe  is  infinite, 
or  they  concern  the  estimation  situation  for  a  large.  In  this  article  we 
concentrate  on  the  practical  cases  of  finite  saaples,  and  smaller  values  of  a. 
Ooodness-of-flt  tests.  Suppose  a  random  sample  x^,x2, . . .  ,xfl  is  given. 

The  goodness -of- fit  tests  are  based  on  the  empirical  distribution 
function  (EDF)  of  the  x^,  and  in  practice  the  statistics  are  calculated  from 
the  values  given  by  the  probability  integral  transform  -  F(x^: 
where  estimates  of  a,  fi  and  a  are  used  in  F(»)  when  the  parameters  are 
not  known.  EDF  tests  have  been  shown  to  be  powerful  in  many  test  situations; 
rival  procedures,  such  as  correlation  tests  or  spaclngs  tests,  often  have 
sero  asymptotic  efficiency  relative  to  EDF  tests  (see,  for  example,  Clbisov 
(1961)  for  ARE  of  spaclngs  tests,  and  McLaren  and  Lockhart  (1987)  for  ARE  of 
correlation  tests). 

The  estimation  techniques  depend  on  the  profile  likelihood  of  the 
data.  They  are  set  out  in  Section  2.  The  tests  of  fit  are  described  in 
Section  3,  and  the  theory  is  given  in  Section  4.  Throughout  the  paper,  the 
plots  of  the  profile  likelihood,  the  estimation  procedures  and  the  tests  of 
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fit  will  be  illustrated  by  rsfsrsnes  to  thrss  data  ssts,  given  in  Appendix  1. 
The  three  sets  ere  as  follows: 

(s)  Data  set  1,  froa  Cox  and  Oakes  (1984,  Table  1.3)  consists  of  10  values  of 

the  nuaber  of  cycles  to  failure  when  springs  are  subjected  to  various  stress 

2 

levels.  For  these  data,  the  stress  level  is  950  N/aa  ,  and  the  values  are  in 
units  of  1000  cycles. 

(b)  Data  set  2  consists  of  15  tines  to  failure  of  air  conditioning  equipaent 
in  aircraft:  these  are  taken  froa  a  table  of  tines  for  several  aircraft, 
given  in  Proachan  (1963,  Table  1),  and  are  the  data  for  aircraft  nuaber  7910. 
This  set  has  been  used  by  Stephens  (1986b)  in  studying  various  tests  for 
exponentiality ,  and  the  conclusion  was  drawn  that  the  tines  con  froa  a 
distribution  with  decreasing  failure  rate  (DFR) ;  if  this  distribution  is 
Veibull,  DFR  iaplies  a  <  1.  Clearly  both  data  sets  1  and  2  have  been 
rounded,  so  they  becoae  discrete,  but  this  nakes  negligible  difference  to  the 
estiaation  procedures  or  the  tests  of  fit. 

(c)  Data  set  3  is  artificially  constructed,  to  Illustrate  the  third  possible 

situation  which  can  occur  (although  aore  rarely)  in  analysing  a  saaple. 

2 

In  order  to  save  space,  we  give  tables  only  for  the  EDF  statistic  A  ,  which 
is  known  in  aany  situations  to  have  good  power.  Tables  for  the  statistics 
and  are  given  by  Lockhart  and  Stephens  (1989). 

Finally  we  reaark  that  tables  for  testing  fit  to  the  larger  faaily  (2) 
above  will  be  published  separately. 


2.  ESTIMATION  PROCEDURES 

2.1  The  different  Cases.  For  the  test  of  #o,  eight  cases  can  be 
distinguished,  according  to  which  paraaeters  in  (1)  oust  be  estiaated;  the 
other  paraaeters  are  assuaed  known.  The  cases  are: 

Case  0:  all  paraaeters  known;  Case  1:  o  unknown;  Case  2:  fi  unknown; 
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Cm*  3:  •  and  fi  unknown;  Casa  4:  a  unknown;  Caaa  5:  a  and  m  unknown; 

Caaa  6:  4  and  a  unknown;  Casa  7:  a,  4  and  a  all  unknown. 

Caaa  Busbars  0,1,2  and  3  correspond  to  the  nuabers  used  for  other 
distributions  Involving  only  location  and  seals  par saa tars ;  see,  for 
axaapla,  Stephans  (1986a).  In  Casa  0,  the  given  by  the  probability 
integral  transform  are ,  on  Hq,  uniform  between  0  and  1;  EOF  tests  for  this 
Caaa  are  given  in  Stephens  (1986a,  Section  4.4).  In  Cases  2,  4,  and  6  above, 
where  a  is  known,  the  transformation  y  -  -  log(x-a)  is  made  and  the 
y- sample  is  tested  to  cosm  from  the  extreme* value  distribution 

F(y)  -  exp  [  -exp ( *  (y-o'  )/4' )  ] ,  -  -  <  7  <  -  (4) 

(Throughout  the  paper,  log  refers  to  natural  logarithm) .  EOF  teats  are  given 
in  Stephens  (1977,  or  1986a,  Section  4.11).  Hare  the  relationship  between  (1) 
and  (4)  is  that  o'  -  -  log4  and  4*  ”  1/a. 

Thus  in  this  article  it  is  necessary  to  give  testa  only  for  Cases  1,  3, 

5  and  7,  where  o  must  be  estimated,  and  possibly  also  4  and  m.  As  was 
stated  in  Section  1,  this  is  usually  done  by  maxima  likelihood,  and  we  first 
examine  the  likelihood  equations. 

Suppose  L(a,fi,a)  is  the  likelihood  fox  a  sample  x^ . xq  from  (1). 

The  log- likelihood  function  is  then 

X  -  log  L( o,4,m)  -  n  log  a  +  (m-1)  Slog(x^-o) 

-  am  log  4  -  SI(xi-o)/4)a  (5) 

where  sums  are  for  i  from  1  to  n.  From  (5)  the  likelihood  equations  are 

-  n/m  +  S  log^-a)  *  (S(x1*a)a  log(x1*a))/4“  - 

(log  4)(«  -  r(xi-o)“/4B)  -  0  (6) 

§£  -  (m-1)  Z(xi*o)'1  *  (a  Z(x1*a)a*1)/4a  -  0  (7) 
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0 


(8) 


Co  give  tho  estimate  of  0  from  Che  estimated  a  end  the  known  or  estimated 
value  of  a. 

The  problems  of  esciaatlon  arias  because ,  when  a  or  a  is 


less  chan  1,  Che  likelihood  can  be  aade  inf  inice  by  seeking  a  -  ,  where 

x,..  <  x.,.  ...  <  x,  .  are  the  order  sCatisCics  of  Che  saaple.  This  is 
clearly  a  biased  esCiaaCe  of  a,  and  we  propose  a  better  estimate  below. 


2.2  Cases  1  and  3;  a  known. 

Suppose  the  known  value  of  a  is  a^;  we  muse  discuss  Che  two  cases  (a) 
aQ  >  1  (b)  0  <  <  1.  (Note  that  when  Bq  -  1,  the  distribution  (1) 

reduces  to  the  exponential  distribution  with  unknown  origin;  tests  for  this 
distribution  have  been  given  by  Spinelll  and  Stephens ,  1986) . 

(a)  Suppose  Bq  >  1.  Then  in  Case  1,  with  0  known,  equation  (7)  is 
solved  for  o.  In  Case  3,  (9)  is  solved  for  a  and  then  (11)  for  0,  using 
a.  It  is  easy  with  computers  to  solve  for  these  estlaates:  a 
straightforward  procedure  is  to  start  with  an  estimated  a(  -  •  c,  where 

e  is  very  small,  and  to  decrease  o^  steadily  until  a  solution  is  found. 
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(b)  Suppose  0  <  «q  <  1.  Ths  M.L.  estimate  of  a  is  now  a  -  x^;  it 
is  «  biased  estimate  and  gives  an  infinite  likelihood.  We  propose  the 
following  alternative  estimates.  Let  k  -  l/m^.  In  Case  1,  take 


For  Case  3,  proceed  iteratively  as  follows:  Suppose  fi( o,m)  is  the 
solution  of  (11),  and  start  with  Oq  -  and  0^  -  0(aQ,ao).  Then  take 

0^  -  ^(or,m^)f  for  r-0,1, .. .until  «r,  0f 
converge  to  estimates  a  and  0.  Many  studies  indicate  diet  this  procedure, 
which  we  call  procedure  A,  converges  rapidly  to  estimates  which  usually  give 
a  better  fit  than  the  usual  a  -  x^,  and  in  Case  3,  0  -  0(a,aQ) .  Procedure 
A  is  illustrated  by  Example  2  below. 

Example  1.  Case  1.  Consider  data  set  1,  in  Table  1.  To  illustrate  Case  1, 
suppose  values  0  -  70  and  m  -  2.1  are  assumed  known.  (These  are 
reasonably  consistent  with  the  values  eet:i«e ted  in  Case  7  below).  The  M.L.K. 
of  a  is  105.31,  from  equation  (7). 

Example  2,  Casa  3.  Consider  data  set  2  in  Table  1,  and  suppose 
a  -  0.5,  say.  Then  usual  estimates  are  a  -  -  12.0,  and  0  -  69.42, 

from  (11).  Procedure  A  converges  rapidly  to  give  o  -  11.69  and  0  -  70.45. 
2,3  5  end  7. 

In  these  cases,  a  and  a  are  both  unknown.  Again  the  likelihood  can 
be  made  infinite,  by  allowing  a  ■  x^j,  and  using  any  estimate  a  such  that 
a  <  1.  Thus  a  cannot  be  made  precise.  We  shall  suggest  estimates  based  on 
the  profile  likelihood.  L*(at) ,  abbreviated  L*.  This  is  the  likelihood  L 
maximised  with  respect  to  0  and  a,  for  a  fixed  a£.  Thus 
L*(ae)  -  L(ae,le,ae)  where,  for  s  given  ac,  for  Case  5,  0  is  the  known 

0,  and  then  a  is  the  solution  of  (6);  for  Case  7,  a  and  0  are  the 

t  fc  C 

solutions  of  (9)  and  (11).  Case  5  (where  only  0  is  known)  occurs  very 


rarely  in  practice,  and  froa  hare  on  w«  shall  diacuaa  estimation  only  for  the 
important  Casa  7,  where  all  throe  parameters  are  unknown. 

Casa  7.  tfe  consider  the  possible  forms  which  a  plot  of  Z(ac)  -  log  L  (at), 
against  a£,  might  take.  Suppose  Z(a£)  is  abbreviated  to  Z.  There  appear 
to  be  only  three  possible  types  of  plot,  illustrated  by  Figures  la,  lb  and 
lc.  These  are  the  plots  for  data  sets  1,  2,  and  3  respectively. 

Figure  la  It  may  be  seen  that  only  Figure  la  has  a  local  maximum,  and  only 
with  this  plot  is  there  a  true  ML  solution,  occurring  for  the  value  of  a 
equal  to  the  a£  at  the  maximum;  then  m  can  be  found  from  (9)  or  (10) ,  and 
0  from  (11)  using  a  and  a.  The  minimum  for  Z  which  occurs  in  Figure  la 
gives  a  saddlepoint  for  the  likelihood.  Figure  la  is  most  likely  to  arise 
when  the  true  m  >  1;  if  the  sample  were  Infinite,  it  would  certainly  be  the 
plot  arising  when  m  >  1.  Rocketts,  Antle  and  Kllmko  (1974)  have  conjectured 
that  Figure  la  occurs  with  probability  converging  to  1  as  n  -»  •  when  m  >  1. 
They  noted  that  the  existence  of  a  local  maximum  implies  the  existence  of  a 
saddlepoint.  Smith  (1985)  shows  inter  alia  that  when  m  >  1,  there  is,  with 
probability  approaching  1  as  »■>•,*  root  of  the  likelihood  equations 
which  is  a  consistent  local  maximum  of  the  likelihood. 

Figure  lb  In  Figure  lb,  Z  steadily  decreases  as  c*c  decreases;  this  figure 
occurs  for  finite  samples  with  increasing  probability  as  m  becomes  smaller; 
for  an  infinite  sample  we  believe  that  it  would  be  certain  to  occur  whenever 
a  <  1.  The  maximum  likelihood  estimate  of  a  would  then  be  giving  an 

infinite  likelihood  for  m  <  1.  This  estimate  is  clearly,  biased,  and  we 
recommend  the  estimate  a  given  by  Procedure  lb  in  Section  2.6  below;  this 
gives  efficient  estimates,  with  a  less  biased  estimate  of  a. 

Figure  lc  Here  there  is  a  minimum  for  Z  (again  this  gives  a  saddlepoint 
for  the  likelihood),  but  there  is  no  maximum  and  therefore  no  ML 
solution.  This  figure  would  effectively  never  occur  if  the  sample  were 
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infinite,  (see  Saith,  1985)  but  with  finite  staples  it  arises  with  increasing 
probability  as  a  increases  (see  Table  4  below).  The  graph  of  Z  tends  to  a 
horizontal  asyaptots  as  a  -»  -  •,  suggesting  that  the  MLE  of  a  -*  in 

C  t 

turn  a  and  0  both  approach  •.  The  Veibull  distribution  then  takes  the 
limiting  fora  discussed  in  the  introduction,  naaely  the  extreae -value 
distribution  (3);  the  relationship  between  the  paraaeters  in  (3)  and  thoae  in 
(1)  is  that  b  -  lia  0/a  and  a  -  lia(a+£),  as  a  -»  0  -*  »,  and  a  -»  •. 

Clearly  if  it  could  be  decided  immediately  from  a  data  set  that  Figure  lc 
would  arise,  there  would  be  no  need  to  try  to  find  the  MLE  for  (1),  and  (3) 
could  be  fitted  directly.  In  Section  2.5  we  show  how  this  My  be  done  by 
calculating  a  quantity  A  whose  sign  iaMdiately  detects  Figure  lc,  but 
first  we  discuss  the  M.L.  equations  (9)  and  (10). 

2.4  Solution  of  the  Maximal  likelihood  equations  (Case  7). 

Maxima  and  ainiM  in  Figures  la,  lb  and  lc  correspond  to  solutions  of  (9) 
and  (10)  for  a  and  a,  and  then  (11)  for  0.  Let  o£  be  a  trial  value  for 
a,  and  let  mg  be  the  solution  of  (9)  and  a^  the  solution  of  (10);  these 
My  be  found  Iteratively.  If  a  plot  of  a^  against  a£  intersects  a  plot 
of  a^Q  against  o^,  we  have  the  solutlon(s)  required.  Plots  of  a^  and  a^Q 
against  a£  are  shown  in  Figures  2a,  2b  and  2c;  they  are  the  plots  for  data 
sets  1,  2  and  3  respectively.  There  are  three  possibilities;  no 
intersection,  so  no  MLE,  as  in  Figure  2b;  one  intersection,  also  no  MLE,  so 
fit  the  extreme -value  distribution,  as  in  Figure  2c;  or  two  intersections , 
with  MLE  corresponding  to  the  solution  with  sMller  o,  as  in  Figure  2a. 

2.5  Detection  of  Figures  lc  and  2c  We  first  give  a  procedure  to  decide  if 
Figures  lc  and  2c  are  appropriate  for  a  given  data  set,  without 

actually  plotting  thea.  In  Figures  2a,  2b,  2c  it  appears  that  the  solutions 
for  Bg  and  tend  to  parallel  lines  as  c»t  -  <*>;  this  is  shown  to  be 


11 


true  in  Appendix  2. 


Suppose  A  is  the  liaiting  "gap"  between  these  lines,  that  is, 

A  “  lin  (a^Q-a.);  the  value  of  A  is  found  as  follows.  Let  x  -  Zx./n,  and 

V‘* 

2  r 

s  -  Hx^/tr,  also  define  Tf  -  Z(x^)  exp(-yx^);  in  these  expressions  the 
suas  run  for  i-l,...n.  Let  7  be  the  solution  of 


The  value  of  7  can  easily  be  found  by  iteration,  starting,  for  exaaple,  with 
7-1  in  the  right-hand  side  terns  Tq  and  T^.  Define 

D  -  iT0  +  7(T2-i  Tx)  ;  (13) 

7  is  the  liaiting  slope  of  the  lines,  and  A  is  given  by 

A  -  (x  Tq  -  7<*T0  -  T2)/2)/D  (14) 

It  is  clear  that  a  negative  value  of  A  laplies  Figure  2c,  and  then  the 
Veibull  fit  should  be  abandoned  in  favour  of  the  extreae-value  fit  (3)  •  note 
that  the  x- values  are  fitted  directly,  and  logarithms  are  not  taken  first,  as 
in  the  test  with  known  a  In  Section  2.  The  detection  of  Figures  lc  and  2c 
by  A  is  equivalent  to  use  of  a  discriminant  L,  proposed  by  Cheng  and  lies 
(1990)  although  they  do  not  discuss  the  behaviour  of  equations  (9)  and  (10) 
in  the  sane  detail. 


Exaaple  3,  Case  7.  Data  set  3  has  been  constructed  to  illustrate  Figure  lc. 
When  the  extreae-value  distribution  (3)  is  fitted  (see,  for  example, 
Stephens,  1977,  or  1986a,  Section  4.10),  the  parasMter  estimates  are 


a  -  1.335,  and  b  -  0.343.  In  Section  3,  we  describe  the  test  of  fit  for 


this  distribution.  It  will  not  be  the  sane  as  that  described  in  Stephens 
(1977;  1986a,  Section  4.10)  because  it  aust  be  aade  conditionally  on  the 
occurrence  of  Figure  lc. 


2.6.  Flgnfi  lb  «nd  2b.  Suppose  A  is  positive,  so  that  either  Figure  Is 
(2s)  or  lb(2b)  is  sppropriste,  and  suppose  s  plot  gives  Figure  lb.  The 
eonventionsl  estimate  will  be  a  -  x^,  with  an  Infinite  likelihood,  but 
again  we  propose  Procedure  A,  adapted  for  unknown  m,  to  give  a  less  biased 
estimate  of  a,  and  also,  almost  always,  a  better  overall  fit. 

Suppose  ap,  0p,  ap  are  estimates  (we  omit  the  A  symbol)  at  iteration 
r:  find  estimates  0^^,  as  follows: 

Ir 

(a)  Let  -  0^X1  ,  where  k  -  l/mr. 

(b)  Then  solve  (6)  for  ,  using  a  - 

(c)  Use  (11)  to  give  0^,  using  and  a^. 

Iteration  of  steps  (a)  to  (c)  continues  until  the  accuracy  required  for 


a  is  obtained.  Initial  estimates  a^  and  0q  may  be  found  by  setting 
°0  "  x(l)  *nd  continuing  with  steps  (b)  and  (c)  above,  but  using  only  the 


n-l  viluu  x(2).  x(3),...x(n). 
a,  0  and  a  for  Figure  lb. 


The  final  estimates  will  be  the  estimates 


Me  A,  Case  7 


Figure  lb  is  the  plot  for  data  set  2;  the  initial 


estimates  a  ,  0  a  are  then  12.0,  101.01  and  a  -  0.795;  with  Procedure  A, 
o  o  o 

7  iterations  give  final  estimates  o  -  9.313,  0  -  93.50,  at  -  0,763.  Here 


iteration  was  stopped  when  two  successive  values  of  a  differed  by  less  than 


0.001.  We  shall  see  in  Section  3  that  the  second  set  of  estimates  gives  a 


better  fit  to  distribution  (1). 


2.7  Figures  la  and  2a:  Example  5,  Case  7:  Finally  we  turn  to  solutions  based 
on  Figures  la  and  2a.  These  can  almost  always  be  found  very  straightforwardly 
by  solving  (9)  and  (10)  by  Iteration  and  searching  for  the  crossing  in 
Figure  2a  with  smaller  a^.  For  data  set  1,  the  maximum  in  Figure  la,  and  the 
crossing  in  Figure  2a,  occur  at  o  -  99.02,  with  a  -  2.38  and  0  -  78.23. 
(The  saddlepoint  is  at  a£  -  116.960,  a  -  1.008). 

2.8  Comment .  When  Figure  lb  occurs,  we  have  found,  from  many  Monte  Carlo 
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studies ,  that  Proesdurs  A  always  converges ,  and  gives  a  better  fit  to  (1); 
then  m  is  snail.  For  larger  values  of  the  true  n.  Procedure  A  also 
converges  for  sons  data  sets,  but  there  is  also  an  MLB  because  the  plot  is 
like  Figure  la.  In  this  case,  the  MLE  should  be  taken. 

If  one  is  seated  before  a  graphics  terminal ,  it  is  not  difficult  to 
plot  the  profile  likelihood  and  decide  between  the  various  Figures:  also  to 
plot  Bg  and  and  find  the  estimates.  It  is  much  sore  difficult  to 

autonste  the  procedure  for  a  computer.  The  value  of  A  in  Section  2.5  will 
decide  whether  or  not  Figures  lc  and  2c  obtain:  the  problem  is  to  decide 
between  Figures  la  (or  2a)  and  lb  (or  2b) .  Ve  suggest  that  be  started  so 
close  to  that  is  almost  1.00  (for  the  examples  given,  this 

nay  mean  -  o£  of  the  order  of  10"*®).  Then  make  a£  smaller  so  that 

ib^q  -  JBg  also  gets  smaller,  and  either  passes  through  sero  (  so  that  a 
saddlepoint  exists  and  we  have  figures  la  and  2a) ,  at  is  clearly  seen  to 
approach  its  limit  A  without  passing  through  sero  (and  we  have  Figures  lb 
and  2b).  In  the  former  case,  once  the  saddlepoint  has  been  found,  the  steps 
in  decreasing  o^  can  be  made  larger  till  the  MLE  is  found  when  is 

zero  for  the  second  time. 


3.  C00PHE8S-0F-FIT  TESTS. 

In  this  section  the  EOF  tests  are  described.  The  null  hypothesis  is 
Hq:  the  random  sample  x^,  x2»”*»x„  comes  from  distribution  (1). 

(a)  Find  the  estimates  of  unknown  parameters  as  described  above,  and  make  the 

transformation,  for  1  -  1,2 . n,  -  F(x^;  using  the  estimates 

where  necessary. 

(b)  Calculate  statistics  A2,  V2 ,  U2  am  follows: 

A2  -  -  n  -  (l/n)Z(2i-l)[log(*(1))  +  log(l  - 
V2  -  2{z(1)  -  (2i-l)/(2n)}2  +  l/(12n) . 
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i  -  1  to  n,  and  log  means  natural 


IF2  -  V2  -  n(z  -  0.5)2,  where 
i  -  £ z^j/n,  and  where  subs  are  for 
logarithms . 

2 

(We  shall  illustrate  using  only  A  in  the  Examples  below). 

(e)  Let  e  -  1/a,  or  1/a  when  a  is  estimated.  Eater  Table  1,  using  the 

subtable  for  the  appropriate  ease.  When  a  or  a  >  2,  we  have  0  <  c  <  0.5 

and  the  subtable  is  entered  on  the  line  corresponding  to  c;  when 

a  or  a  is  less  than  or  equal  to  2,  so  that  e  t  0.5,  the  last  line  of  the 

subtable,  labelled  c  -  0.5,  should  be  entered.  Bq  is  rejected  at 

significance  level  p  if  the  statistic  used  is  greater  than  the  value  given 

for  level  p.  The  table  has  been  given  using  c  rather  than  a  because 

linear  Interpolation  for  c  will  give  good  accuracy.  In  all  the  tables,  the 

given  points  are  for  the  asynptotic  distributions  of  the  statistics:  however, 

they  can  be  used  with  good  accuracy  for  saallar  values  of  n  (say  n  2  10) ; 

for  n  <  10  a  goodness -of- fit  test  would  in  any  case  have  very  little  power. 

Example  1  Case  1  (continued)  For  data  set  1,  with  a  assumed  to  be  2.1,  (so 

c  -  1/a  -  0.476).  fi  -  70,  and  a  -  105.31,  the  value  of  A2  is  0.307. 

Reference  to  Table  1  for  Case  1,  with  entry  at  c  •  0.476,  shows  no 

significance  at  the  50%  level,  so  the  Welbull  fit  is  good. 

Example  2.  Case  3  (continued)  For  data  set  2,  with  a  -  0.5,  a  «■  11.69  and 

fi  -  70.44  as  in  Example  2,  we  have  A  -  .754.  Since  c  -  1/a  -  2.0,  greater 

2 

than  0.5,  the  table  for  Case  3  is  entered  on  the  last  line  (c  -  0.5);  A  is 
not  significant  at  the  20%  level,  so  the  fit  is  good. 

Example  3.  Case  7  For  data  set  3,  the  extreme -value  distribution  (3)  must 

be  fitted  and  estimates  are  a  -  1.335  and  h  -  0.343.  After  the  Probability 

Integral  Transform  z ^  -  F(x^),  where  F(x)  is  given  by  (3),  the  value 
2 

of  A  is  0.561.  This  value  should  then  be  referred  to  Table  1,  Case  7, 
with  c  -  0,  corresponding  to  a  -  •.  The  p -value  is  approximately  0.08. 


Co— at.  This  tsst  should  not  bs  confused  with  tho  usual  tost  for  tbs 

i 

extras- value  distribution  with  unknown  psrsastors,  givon  by  Stephens  (1977; 

1986s,  Section  4.10).  The  test  given  in  the  previous  section  is  asde 

conditionally  on  the  occurrence  of  Figure  lc  for  the  profile  likelihood  plot, 

end  Table  3  given  here,  with  c  -  0,  is  then  the  relevant  table.  The  tests 

given  by  Stephens  ere  those  for  the  situation  where  it  is  intended  from  the 

start  to  fit  ea  extreae -value  distribution  to  the  date  set. 

Example  4  Case  7  (continued)  For  data  set  2,  assuming  all  three  parameters 

estimated  in  the  conventional  way,  giving  a  -  12.00,  £  -  101.01,  s  -  0.795, 

2 

we  have  A  -  0.703.  When  Procedure  A  is  used,  the  estimates,  given  in 
Section  2.10,  are  now  a  -  9.313,  fi  -  93.50,  a  -  0.763,  and  the  test  statistic 
is  J?  -  0.54,  indicating  a  batter  fit.  Since  c  *  1/a  -  1.31  is  greater  than 
0.5,  the  Table  for  Case  7  is  entered  on  the  last  line  (c  -  0.5).  The 
significance  level  is  then  0.17. 

Example  3  Case  7  (continued)  For  data  set  1,  with  all  three  parameters 
estimated  as  in  Example  3  above,  namely  a  -  99.02,  £  -  78.23,  and  m  -  2.38, 

-  0.260.  Table  1  is  entered  at  c  -  1/a  -  0.420.  The  above 
test  value  is  not  significant  at  the  50%  level,  so  the  fit  is  very  good. 

4.  ASYMPTOTIC  THEORY  OF  EDF  TESTS. 

4.1  Asymptotic  distributions 

In  this  section  the  asymptotic  theory  of  EDF  tests  is  summarized.  The 

calculation  of  asyaptotic  distributions  of  EDF  statistics  follows  a 

well-known  procedure  (see,  for  example,  Durbin,  1973;  Stephens,  1976).  It  is 

based  on  the  fact  that  y_(z)  -  ^n{F  (z)  -  x) ,  where  F(x)  ia  the  EDF  of 

n  xi  n 

the  z-set,  tends  to  a  Gaussian  process  y(z)  as  n  •»  «,  and  the  statistics 
are  functionals  of  this  process.  The  mean  of  y(z)  is  zero:  we  need  the 
covariance  function  p(s,t)  ■  E{jr(s)  y(t)J.  When  all  paraaeters  are  known 
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(Cm*  0),  this  covariance  la  -  ain(jr,  e)  -  at.  Whan  par  an*  tars  are 

estimated,  die  covariance  will  depend  on  those  which  are  estimated;  if  the 
method  of  estimation  is  efficient  the  covariance  will  not  depend  on  true 
values  of  location  or  scale  par  asm  tars  a  and  fi ,  hut  will  depend  on  the 
shape  parameter  a.  We  Illustrate  the  calculation  for  Cm*  7,  the  most 
difficult  CMC. 

Suppose  the  parameters  are  components  of  a  vector  I:  -  a,  9^  m  4. 

#3  -  a,  and  (Case  7)  suppose  all  three  are  unknown.  Let  F(x;#)  now  denote 
the  distribution  and  let  f(x;l)  be  the  corresponding  density. 

Suppose  a  vector  g(s),  with  components  g^(*)  is  constructed  as  follows: 

gl(t)  "  '  1“1’2*3  •  <l5> 

where  the  right  hand  side  is  written  as  a  function  of  s  using  the 
transformation  s  -  F(x;#).  Let  (f(a)}#  denote  the  transpose  of  g(s). 

Let  D  be  the  (symmetric)  matrix  with  entries 

SLj  -  £(-a2logf(x;#)/S#1S#^),  i.J  -  1,2,3,  where  E  denotes  expectation,  and 
let  £  be  the  inverse  of  D.  Then,  for  Cmo  7, 


P7(s,e)  -  po(s,t)  -  (g(s) )'£*(*)  .  (16) 
From  (23)  above,  the  components  of  g(s)  became,  after  some  algebra,  and 
using  F  for  F(x;0): 


*2<*>  * 


*3<*>  - 


SF  _  _  (1-s) 

dm  “  m 


(log(l-s))  log( ~log(l~s) ) 


Also,  for  Cue  7,  and  for  m  >  2,  D  has  the  top  right  terms: 


When  £  is  calculstsd  and  |(<)  and  S  ara  inaartad  into  (16)  ^(a,C) 
will  ba  indapandant  of  a  and  fi.  Whan  ad  2,  Cha  M.L.  estimate  a  of  a 
is  supor  officiant  in  tha  a ansa  of  Darling  (1955)  and  than  tha  covariance  will 
not  noad  tha  first  tarn  in  g(s)  and  corraaponding  tans  in  D;  tha 
asyeptotics  ara  the  sane  as  if  a  wars  known,  that  is,  for  Casa  6.  Thus 
f(«)  -  (fjU))  fj(«))  and  D  is  as  in  (18)  but  with  tha  top  row  and  first 
colusn  mowed. 

For  other  cases,  ^(«,t)  for  Casa  k  is  calculated  using  only  those 

components  in  ^(i)  which  correspond  to  unknown  1^,  with  tha  corresponding 

entries  in  D,  before  this  is  inverted  to  give  the  Z  used  in  (16).  The 

Craa6r-won  Miaes  statistic  ^  is  based  directly  on  tha  process  y(x) ,  while 

d2  is  based  on  tha  process  a(s)  -  y(*)/(s(l-*))x/*j  asyeptotically  W2  and 
2 

A  are  given  by 

1  1 

V*  -  J  y2(*)d*  and  A2  -  |  a2(s)d*  . 

0  0 

Tha  asyuptotic  distributions  of  both  statistics  are  a  sun  of  weighted 
2 

independent  Xj  variables;  tits  weights  aust  be  found  froa  the  eigenvalues  of 

an  integral  equation  with,  for  V2,  the  pk(s, t)  for  case  k,  as  kernel. 

2 

For  A  one  aust  find  the  ^(>,t)  of  tha  a(z)  process.  Once  tha 
weights  are  known,  tha  percentage  points  of  the  distributions  can  be 
calculated  by  Xahof's  aethod.  The  techniques  are  straightforward  once  the 
V*  , t)  are  known,  rad  we  onit  die  details:  they  are  given  in  Lockhart  and 
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SCifhini  (1919) . 


5.  VUITHOL  WttlH 

5.1  Aiyptotlc  dlmlktiflM  and  Moot*  Carlo  itndUt  The  results  of  Smith 

(1985)  esn  bo  used  to  establish  rigorously  that  has  the  asymptotic 

distributions  calculated  as  above,  for  any  estimator  found  by  a  method  of 

estimation  having  the  properties  of  Theorem  3  of  Smith.  In  an  unpublished 

note  the  present  authors  have  demonstrated  the  existence  of  such  an  estimator 

by  establishing  that  it  is  possible  to  select  the  "correct"  local  maximum  of 

the  likelihood  in  the  event  that  there  are  several.  However,  we  emphasize 

that  op  Monte  Carlo  data  set  haa  aver  been  found  in  which  two  such  local 

2 

maxima  exist.  MS  should  also  note  hare  that,  for  A  ,  we  are  unable  to  give  a 
rigorous  derivation  for  the  asymptotic  distribution;  this  problem  arises 
whenever  parameters  are  estimated  in  goodness -of -fit  tests  of  the  Anderson- 
Darling  type.  The  asymptotic  distribution  is,  however,  confirmed  by  extensive 
Monte  Carlo  studies.  These  Monte  Carlo  studies  were  undertaken  to  confirm 
various  features  of  the  above  estimation  and  tasting  procedures.  Typically, 
they  involved  10000  samples  of  sixes  n  -  10,  20,  30,  40,  60,  100,  snd  200,  and 
with  values  of  a  from  0.1  to  10.  The  studies  firstly  confirmed  the 
conjecture,  proposed  by  other  authors  also,  that  plots  of  m^  and  m^  will 
not  cross  more  than  twice,  so  that  a  local  maximum  of  the  likelihood  occurs, 
if  at  all,  only  once.  This  is  clearly  important  in  knowing  what  to  search  for 
when  the  parameters  are  to  be  estimated.  The  studies  also  confirmed  the 
success  of  Procedure  A  in  providing  estimates  to  give  a  better  fit,  noted 
earlier. 

5.2.  Frequency  of  Figures  la  lb  and  lc.  These  studies  also  revealed  how  the 
relative  frequency  of  Figures  la,  lb  and  lc  will  depend  on  n  and  on  m.  As 
n  grows  larger,  Figure  lc  becomes  less  and  less  likely  and  as  m  becomes 


smaller,  for  fixed  n,  the  relative  frequency  changu  froa  Figure  la  towards 
Figure  lb.  Thasa  results  are  Illustrated  in  Table  4. 

The  results  of  Saith  (1985)  guarantee  that  Figure  lc  arises  only  with 
probability  tending  to  0  as  n  tends  to  infinity  for  any  fixed  finite  value 
of  a  in  (1).  If  the  data  are  saapled  froa  the  two  paraaeter  extraaa  value 
distribution,  however,  (that  is,  (1)  with  s  4  •)  Figure  le  nay  be  expected  to 
arise  about  half  the  ties.  This  can  be  explained  as  follows.  The  saapllng 
can  be  regarded  as  froa  the  family  (2),  with  e  -  0,  and  sample*  will  give 
sonatinas  c  >  0  (Figure  la),  and  sonatinas  c  <  0  (Figure  lc).  This  problen 
of  unobtainable  estinates  within  the  desired  f sadly  also  arises  in  testing  for 
the  von  Mises  distribution  on  the  circle  with  known  direction  of 
concentration.  In  Lockhart  and  Stephans  (1985)  it  was  suggested  that  there 
too  an  expansion  of  the  nodel  could  overcoaa  this  problen.  Slailarly  here, 
tasting  for  the  Generalized  Extraaa  Value  distribution  (2)  rather  than  the  3 
paraaeter  tfeibull  distribution  (1)  will  aliainate  the  awkwardness.  Tests  for 
the  enlarged  nodal  (2)  will  be  presented  in  a  future  paper. 

5.3  Convergence  of  distributions  of  IDF  statistics  For  finite  n,  the  Monte 
Carlo  studies  show  that  the  distributions  of  W*,  U*  and  converge  rapidly 
to  the  asymptotic .  This  is  similar  to  the  behaviour  of  these  statistics  in 
other  test  situations  and  tables  will  not  be  given.  The  asymptotic  points  in 
Table  1  can  be  used  with  good  accuracy  for  n  k  10. 

EDF  statistics  are  known  to  provide  powerful  tests  for  many 

distributions;  the  powers  naturally  depend  on  the  alternatives  considered, 

and  a  study  is  being  made  on  power  properties  for  the  various  alternatives  to 

the  Weibull  usually  encountered.  On  the  whole,  with  the  limited  power 

2 

results  at  present  available,  the  statistic  A  above  is  suggested  as  the 
preferred  statistic  for  overall  Veibull  testing. 
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APfEHDIX  1:  BATA  SETS 


S«t  3.  Artificial  data. 

0.273  0.468  0.504  0.535  0.617  0.804  0.932  1.034  1.289  1.293 

1.294  1.376  1.399  1.407  1.422  1.497  1.521  1.542  1.685  1.737 
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APPENDIX  2 


(«)  Asyaptotas  for  a^  and  a^  (Sactlon  2.5) 

Supposa,  as  a  th«  ralacion  a  -  ya  +  +  fc^/o,  vhara  7, 

and  and  b^  ara  Co  ba  datarainad  balow  for  aaeh  of  a^  and  a^g. 

Than  va  hava 


1  - 


r)  -«p|-  + 

■U 


(Al) 


Vxi  t^3  +  +  Vii 

—  + - 5 - 


- y*,  -y*i  9 

a  Aaxp(-Z),  aay,  -  a  x  11  -  IZ)  +  [2*/ 2)] 


-7*i 

■ 

,  ** i /2  +  b0*i  . 

3  2 

fr****J 
- £ - - - 

2  4  .22 

a 

a 

V*I 


n. 


Hanca 

va  gat 
n 


_  r  ’JX1 

•acting  Tr  -  2.  xi  a  ,  vhara  suaa  ara  for  1  froa  1  Co 


*3 

B, 


x  (1  -  x^o)*  -  r0  +  +  -4  vhara 

I-l  a 

d2  -  -  [tT2/2  bQTx]  and 

B2  "  *  <?V3  +  boV2  +  5iri)  +  “^V8  +  W2  +  b0  V2  ’ 

Also  Cj  -  X(l-x1/o)B'1  -  Tq  +  ~  +  Bj/o2  ,  vhara 

-  [7T2/2  +  (Bq-D^]  and 

-  (7  r3/3+(b0-i)  r2/2  +  5^)  +  72r4/8  +  (b0-i)7r3/2  +  (b2-2b0+i)r2/2 


V*  now  use  tho  above  results  in  aquation  (10).  This  nay  ba  writtan 
Z(1  •  xL/a)'1  -  n  Z(1  -  xL/a)a'1  . 


Carrying  thraa  tarns  in  aach  expression,  wa  have 


■». +  V  1  * 


+  -J  -  [t*  +  bQ  + 

a  ' 


bl)  f  *3 
—  rn  +  ~ 

°Jl°  •  «■ 


where  x  -  Zx^/n  and  s  -  Sx^/n.  The  coaff icianta  of  the  constant  tarn  give 
V*0  +  <b0*1)20  +  TAj  -  7^3  +  b0r0;  than  7(irQ  ♦  \  -  V  -  rQ.  But 
dj  -  Aj  -  -  T^,  so  we  have  for  7, 

1  -  (i  -  ^/Tg)  .  <A2) 

Next,  equating  coefficients  of  1/a,  we  have 

7s  rQ  +  7  i  ^  +  7*2  ♦  <V1)iT0  +  (b0'1)A2  *  blT0  “  7  *3  +  Vs  +  blT0; 

then  7(5  rQ  +  xAj  +  B2  -  B3)  +  ^  ’  *3*  “  **0  *  *2  “  0  *  (A3) 


b2- 


b3- 


72r.  bn7T-  b«T- 


[*3  V2  .  .  _  )  .  *  r4  b0*3  b7l 

■  [-y  +  t—  +  6irlJ  +  —  +  ~~2~  +  “T 


P-2 


0-1)T2  .  .  .  1  .  l\  .  .  »S  ’  2Vl>T2 

-X - +  b,T,  +  — a-  +  - s -  + - x - 


Bj  -  Bj  -  -j-  +  b^Tg  -  r2  .  Recall  also  that  Aj  -  than. 


substitution  in  (A3)  gives 


b0(ir0  +  7  T2  -  t^J  +  2_  (r3-  ir2)  +  j(stq  -  r2/2) 


-  xro-0 


Solution  of  (A2)  for  y,  and  (A4)  for  bg,  now  callad  gives  tho 

coofficionts  in  tho  aayaptoto  m^Q  -  ya  +  f  glvon  in  Soction  2.5. 

Vo  now  eonoldor  equation  (9),  which  nay  bo  written 

S(1  -  x^o)*  -  «|r  logd-x^o) (1  -  x^o)-  -  Z  log(l.x1/a)js(l-x1/o)“| 
Expanding,  wo  have 


-  7(^0-^)  +  jbQtxTQ.^)  +  yisTf/ 2  +  a2 i  +  (7r3-r2)/2  +  b0r2]|/o  . 

Equating  tho  constant  torn  gives  Tq  -  yCxT^-T^) ,  tho  sane  as  for  equation 
10.  Thus  tho  asymptotes  for  and  ^  will  bo  parallel.  Further,  the 
coefficient  of  1/a  gives 

2 

7*Tq  _  yTy  j  T- 

Aj  -  bQ(xT0  -  rx)  +  7  -j~  +  7d2x  -  -5-  +  — y~  +  7&or2  • 

■ore  algebra  gives 

-  -  2  f^3  ^2*1  ^0 

b0(xT0  7r2  -  T^x)  +  /jji  -  -|-J  +  -j-2  -  0  (A5) 

Solution  of  (A5)  for  gives  the  constant  6g  in  the  asymptote 

a9  -  ya  +  6g  quoted  in  Section  2.5. 

Difference  A  Finally,  we  have  A  -  S ^  -  Sg  given  by 

(«10-fi9)(ir0  +  7T2-  7xr^)  -  ir0  -  7<sr0-r2)/2  , 

as  in  equation  (14)  of  Section  2.5. 


APPENDIX  3 


In  Section  3.1  it  vu  pointed  out  that  the  aeons  of  the  esyaptotlc 

distributions  of  W2  and  U2  can  be  found  analytically.  For  ccwpleteness , 

we  list  below  six  integrals  which  arise  in  these  calculations. 

1 

Ix  -  J  (l-*)2  {»log(l-s))2<B'1)/*  ds  -  r<3-2/a)/3(3*2/B)  ; 

0 

1 

12  -  J  (l-s)2{log(l-s))2ds  -  2/27 

0 

1 

13  -  J  (l-s)2{log(l-s) )2{log{ -log(l-s) } ]2  ds 

0 

-  2(*2/6  +  y2-  3y  +  2  +  (2y3)log  3  +  log23>/27 

-  0.105618 

1 

14  -  J  <l-s)2{-log<l-s))2‘1/B  ds  -  r<b)/3b 

0 

1 

15  -  J  (l-s)2{-log(l-s)J2‘1/B  log{ -log(l-s) )ds 

0 

-  tr'(b)  -  iog3  r<b)j/3b  . 

where,  in  1^  and  Ij,  b  -  3-1/a; 

1 

Ig  -  J  (l-s)2{log(l-s) J2  log{ -log(l-s) )ds 
0 

-  tr'(3)  -  2  log3}/27 

-  {3  -  2y  -  2  log3)/27  -  -  0.01302 


TABLE  I 

CRITICAL  POINTS  FOR  WEI BULL  CASE  1  FOR  A**2 


c 

0.500 

0.750 

0.850 

0.900 

0.950 

0.975 

0.990 

o.o 

! 0.496 

0.735 

0.915 

1.061 

1.321 

1.590 

1.958 

0.05 

0.484 

0.713 

0.882 

1.019 

1.260 

1.510 

1.851 

0.10 

0.474 

0.693 

0.853 

0.982 

1.208 

1.440 

1.756 

0.15 

0.467 

0.677 

0.830 

0.953 

1.166 

1.385 

1.680 

0.20 

0.464 

0.670 

0.818 

0.936 

1.141 

1.350 

1.631 

0.25 

0.468 

0.674 

0.821 

0.938 

1.139 

1.343 

1.615 

0.30 

0.485 

0.696 

0.847 

0.966 

1.171 

1.377 

1.651 

0.35 

0.517 

0.747 

0.910 

1.039 

1.259 

1.480 

1.774 

0.40 

0.572 

0.840 

1.032 

1.185 

1.448 

1.713 

2.068 

0.45 

0.654 

0.995 

1.249 

1.455 

1.816 

2.189 

2.694 

0.50 

0.774 

1.248 

1.621 

1.933 

2.492 

3.077 

3.878 

CRITICAL  POINTS  FOR  WEI  BULL  CASE  3 

FOR  A**2 

C 

0.500 

0.750 

0.850 

0.900 

0.950 

0.975 

0.990 

•  V, 

0.0 

0.342 

0.472 

0.563 

0.635 

0.757 

0.879 

1.043 

0.05 

0.345 

0.476 

0.568 

0.640 

0.763 

0.886 

1.051 

0.10 

0.349 

0.482 

0.575 

0.648 

0.773 

0.898 

1.065 

0.15 

0.354 

0.491 

0.586 

0.661 

0.789 

0.917 

1.088 

0.20 

0.368 

0.503 

0.602 

0.679 

0.812 

0.945 

1.122 

0.25 

0.374 

0.520 

0.624 

0.705 

0.844 

0.984 

1.171 

0.30 

0.388 

0.544 

0.654 

0.740 

0.889 

1.039 

1.240 

0.35 

0.407 

0.574 

.0.6.94 

0.788 

0.951 

1.116 

1.338 

0.40 

0.430 

0.614 

0.747 

0.853 

1.037 

1.224 

1.478 

0.45 

0.459 

0.667 

0.819 

0.941 

1.156 

1.376 

1.675 

0.50 

0.496 

0.735 

0.915 

1.061 

1.321 

1.590 

1.958 

i 


0.995 


2.243 

2.115 

2.001 

1.909 

1.847 

1.825 

1.860 

1.997 

2.340 

3.086 

4.383 


0.995 


1.167 

1.176 

1.193 

1.219 

1.258 

1.314 

1.394 

1.509 

1.673 

1.909 

2.243 
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1* 
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1 

J  CRITICAL  POINTS  FOR  WEI BOLL  CASE  5  FOR  A**2 

> 


c  J 

0.500 

0.750 

0.850 

0.900 

0.950 

0.975 

0.990 

0.995 

1 

0.0 

.  0.342 

0.472 

0.563 

0.635 

0.757 

0.879 

1.043 

1.167 

0.05  f  0.343 

0.473 

0.564 

0.636 

0.758 

0.881 

1.045 

1.170 

0.10  ; 

0.344 

0.476 

0.567 

0.640 

0.763 

0.887 

1.052 

1.178 

0.15  , 

0.348 

0.480 

0.574 

0.647 

0.772 

0.898 

1.067 

1.195 

0.20  ! 

0.353 

0.489 

0.585 

0.660 

0.789 

0.919 

1.092 

1.225 

0.25 

0.363 

0.505 

0.604 

0.688 

0.818 

0.955 

1.138 

1.278 

0.30 

0.381 

0.532 

0.640 

0.725 

0.872 

1.021 

1.222 

1.377 

0.35  ' 

0.411 

0.583 

0.706 

0.805 

0.977 

1.154 

1.394 

1.580 

0.40 

0.460 

0.673 

0.832 

0.962 

1.192 

1.432 

1.762 

2.017 

0.45 

0.535 

0.828 

1.061 

1.257 

1.611 

1.982 

2.490 

2.883 

0.50 

0.634 

1.059 

1.418 

1.724 

2.276 

2.853 

3.639 

4.247 

C 

CRITICAL  POINTS  FOR  WEI BOLL 

CASE  7 

FOR  A**2 

0.995 

0.500 

0.750 

0.850 

0.900 

0.950 

0.975 

0.990 

0.0  . 

0.292 

0.395 

0.467 

0.522- 

0.617 

0.711 

0.836 

0.931 

0.05 

0.295 

0.399 

0.471 

0.527 

0.623 

0.719 

0.845 

0.941 

0.10 

0.298 

0.403 

0.476 

0.534 

0.631 

0.728 

0.856 

0.954 

0.15  * 

0.301 

0.408 

0.483 

0.541 

0.640 

0.738 

0.869 

0.969 

0.20 

0.305 

0.414 

0.490 

0.549 

0.650 

0.751 

0.885 

D.986 

0.25 

0.309 

0.421 

0.498 

0.559 

0.662 

0.765 

0.902 

1.007 

0.30 

0.314 

0.429 

0.508 

0.570 

0.676 

0.782 

0.923 

1 .030 

0.35  * 

0.320 

0.438 

0.519 

0.583 

0.692 

0.802 

0.947 

1.057 

0.40 

0.327 

0.448 

0.532 

0.598 

0.711 

0.824 

0.974 

1.089 

0.45 

0.334 

0.469 

0.547 

0.615 

0.732 

0.850 

1.006 

1.125 

0.50 

0.342 

0.472 

0.563 

0.636 

0.757 

0.879 

1.043 

1.167 

I 
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TABLE  2 


CRITICAL  POINTS  FOR  WEI BULL  CASE  1  FOR  W**2 


C  0.500  0.750  0.850  0.900  0.950  0.975  0.990  0.995 


0.0 

0.074 

0.116 

0.148 

0.174 

0.222 

0.271 

0.338 

0.390 

0.05 

0.071 

0.111 

0.140 

0.164 

0.207 

0.251 

0.312 

0.358 

0.10 

0.069 

0.106 

0.134 

0.156 

0.195 

0.235 

0.289 

0.332 

0.15 

0.068 

0.103 

0.129 

0.149 

0.186 

0.222 

0.272 

0.310 

0.20 

0.067 

0.101 

0.126 

0.146 

0.180 

0.215 

0.261 

0.296 

0.25 

0.068 

0.102 

0.127 

0.147 

0.180 

0.214 

0.259 

0.293 

0.30 

0.071 

0.107 

0.133 

0.154 

0.189 

0.224 

0.271 

0.306 

0.35 

0.077 

0.117 

0.147 

0.171 

0.211 

0.252 

0.306 

0.347 

0.40 

0.086 

0.136 

0.172 

0.202 

0.254 

0.308 

0.381 

0.436 

0.45 

0.100 

0.164 

0.215 

0.257 

0.332 

0.410 

0.518 

0.601 

0.50 

0.119 

0.209 

0.284 

0.347 

0.461 

0.581 

0.744 

0.869 

C 

CRITICAL  POINTS  FOR 

WEI BULL 

CASE  3 

FOR  W**2 

0.995 

0.500 

0.750 

• 

0.850 

0.900 

0.950 

0.975 

0.990 

0.0 

0.050 

0.073 

0.089 

0.102 

0.124 

0.146 

0.175 

0.197 

0.05 

0.051 

0.074 

0.091 

0.104 

0.126 

0.149 

0.179 

0.202 

0.10 

0.052 

0.076 

0.093 

0.107 

0.130 

0.153 

0.184 

0.208 

0.15 

0.054 

0.078 

0.096 

0.1 10 

0.134 

0.159 

0.191 

0.216 

0.20 

0.055 

0.081 

0.100 

0.115 

0.140 

0.166 

0.201 

0.227 

0.25 

0.057 

0.085 

0.105 

0.121 

0. 148 

0.176 

0.213 

0.241 

0.30 

0.060 

0.089 

0.110 

0.128 

0.157 

0.187 

0.228 

0.260 

0.35 

0. 063 

0.094 

0.117 

0.136 

0.169 

0.202 

0.248 

0.283 

0.40 

0 , 066 

0.100 

0.126 

0.147 

0.183 

0.221 

0.272 

0.311 

0.45 

3.070 

0.107 

0.136 

0.159 

0.200 

0.243 

0.301 

0.346 

0.50 

0.074 

0.116 

0.148 

0.174 

0.222 

0.271 

0.338 

0.390 

30 


CRITICAL  POINTS  FOR  WEI BOLL  CASE  5  FOR  W**2 


C  0.500  0.750  0.850  0.900  0.950  0.975  0.990  0.995 


0.0 

0.050 

0.073 

0.089 

0.102 

0.124 

0.146 

0.175 

0.197 

0.05 

0.050 

0.073 

0.089 

0.102 

0.124 

0.146 

0.176 

0.198 

0.10 

0.051 

0.074 

0.090 

0.103 

0.126 

0.148 

0.178 

0.200 

0.15 

0.052 

0.075 

0.092 

0.105 

0.128 

0.151 

0.182 

0.206 

0.20 

0.053 

0.078 

0.095 

0.109 

0.133 

0.157 

0.189 

0.214 

0.25 

0.055 

0.081 

0.100 

0.115 

0.141 

0.167 

0.202 

0.229 

0.30 

0.059 

0.088 

0.109 

0.126 

0.155 

0.185 

0.226 

0.257 

0.35 

0.065 

0.098 

0.124 

0.144 

0.180 

0.218 

0.270 

0.310 

0.40 

0.073 

0.116 

0.150 

0.178 

0.228 

0.281 

0.353 

0.409 

0.45 

0.086 

0.145 

0.194 

0.236 

0.311 

0.390 

0.498 

0.581 

0.50 

0.102 

0.186 

0.258 

0.320 

0.431 

0.547 

0.706 

0.827 

C 

CRITICAL  POINTS  FOR  WEI  BOLL 

CASE  7 

FOR  W**2 

0.995 

0.500 

0.750 

0.850 

0.900 

0.950 

0.975 

0.990 

0.0 

0.044 

0.062 

0.075 

0.085 

0.103 

0.120 

0.144 

0.162 

0.05 

0.044 

0.063 

0.076 

0.086 

0.104 

0.122 

0.145 

0.163 

0.10 

0.044 

0.063 

0.077 

0.087 

0.105 

0.123 

0.147 

0.165 

0.15 

0.045 

0.064 

0.077 

0.088 

0.106 

0.125 

0.149 

0.168 

0.20 

0.045 

0.065 

0.079 

0.089 

0.108 

0.127 

0.152 

0.170 

0.25 

0.046 

0.066 

0.080 

0.091 

0.110 

0.129 

0.154 

0.174 

0.30 

0.047 

0.067 

0.081 

0.093 

0.112 

0.132 

0.157 

0.177 

0.35 

0.047 

0.068 

0.083 

0.094 

0.114 

0.134 

0.161 

0.181 

0.40 

0.048 

0.069 

0.085 

0.097 

0.117 

0.138 

0.165 

0.186 

0.45 

0.049 

0.071 

0.087 

0.099 

0.120 

0.141 

0.170 

0.191 

0.50 

0.050 

0.073 

0.089 

0.102 

0.124 

0.146 

0.175 

0.197 
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TABLE  3 


CRITICAL  POINTS  FOR  WEI BULL  CASE  1  FOR  U**2 


C 

0.500 

0.750 

0.850 

0.900 

0.950 

0.975 

0.990 

0.995 

0.0 

0.060 

0.090 

0.111 

0.129 

0.159 

0.189 

0.230 

0.261 

0.05 

0.060 

0.090 

0.111 

0.129 

0.159 

0.189 

0.230 

0.261 

0.10 

0.060 

0.090 

0.112 

0.129 

0.159 

0.189 

0.230 

0.262 

0.15 

0.060 

0.090 

0.112 

0.130 

0.160 

0.190 

0.231 

0.262 

0.20 

0.061 

0.091 

0.113 

0.131 

0.161 

0/191 

0.232 

0.264 

0.25 

0.062 

0.092 

0.114 

0.132 

0.163 

0.193 

0.235 

0.266 

0.30 

0.063 

0.094 

0.116 

0.134 

0.165 

0.196 

0.238 

0.270 

0.35 

0.064 

0.096 

0.119 

0.138 

0.169 

0.201 

0.243 

0.276 

0.40 

0.066 

0.100 

0.124 

0.143 

0.176 

0.209 

0.253 

0.286 

0.45 

0.069 

0.105 

0.131 

0.152 

0.187 

0.222 

0.268 

0.304 

0.50 

0.069 

0.105 

0.131 

0.152 

0.187 

0.222 

0.268 

0.304 

CRITICAL  POINTS  FOR  WEI BULL  CASE  3  FOR  U**2 


C 

0.500 

0.750 

0.850 

0.900 

0.950 

0.975 

0.990 

0.995 

0.0 

0.048 

0.070 

0.085 

0.097 

0.118 

0.139 

0.167 

0.187 

0.05 

0.049 

0.071 

0.086 

0.098 

0.119 

0.140 

0.168 

0.189 

0.10 

0.049 

0.071 

0.087 

0.100 

0.121 

0.142 

0.170 

0.192 

0.15 

0.050 

0.073 

0.089 

0.101 

0.123 

0.145 

0.174 

0.195 

0.20 

0.051 

0.074 

0.091 

0.104 

0.126 

0.148 

0.178 

0.200 

0.25 

0.052 

0.076 

0.093 

0.107 

0.130 

0.153 

0.184 

0.207 

0.30 

0.054 

0.079 

0.096 

0.111 

0.135 

0.159 

0.191 

0.215 

0.35 

0.055 

0.082 

0.100 

0.115 

0.141 

0.167 

0.201 

0.227 

0.40 

0.058 

0.085 

0.105 

0.121 

0.149 

0.176 

0.214 

0.242 

0.45 

0.060 

0.090 

0.111 

0.129 

0.159 

0.189 

0.230 

0.261 

0.50 

0.060 

0.090 

0.111 

0.129 

0.159 

0.189 

0.230 

0.261 
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CRITICAL  POINTS  FOR  WEI BULL  CASE  5  FOR  U**2 


c 

0.500 

0.750 

0.850 

0.0 

0.048 

0.070 

0.085 

0.05 

0.048 

0.070 

0.085 

0.10 

0.049 

0.070 

0.086 

0.15 

0.049 

0.071 

0.087 

0.20 

0.050 

0.072 

0.088 

0.25 

0.051 

0.074 

0.090 

0.30 

0.052 

0.076 

0.093 

0.35 

0.053 

0.078 

0.096 

0.40 

0.055 

0.082 

0.101 

0.45 

0.057 

0.085 

0.106 

0.50 

0.057 

0.085 

0.106 

0.900 

0.950 

0.975 

0.990 

0.097 

0.118 

0.138 

0.166 

0.097 

0.1 18 

0.139 

0.167 

0.098 

0.119 

0.140 

0.168 

0.099 

0.121 

0.142 

0.170 

0.101 

0.123 

0.145 

0.174 

0.103 

0.126 

0.148 

0.178 

0.106 

0.130 

0.153 

0.185 

0.110 

0.135 

0.160 

0.194 

0.116 

0.143 

0.170 

0.207 

0.123 

0.152 

0.181 

0.222 

0.123 

0.152 

0.181 

0.222 

CRITICAL  POINTS  FOR  WEI BULL  CASE  7  FOR  0**2 


C 

0.500 

0.750 

0.850 

0.0 

0.043 

0.061 

0.074 

0.05 

0.043 

0.062 

0.075 

0.10 

0.044 

0.062 

0.076 

0.15 

0.044 

0.063 

0.077 

0.20 

0.045 

0.064 

0.077 

0.25 

0.045 

0.065 

0.078 

0.30 

0.046 

0.065 

0.080 

0.35 

0.046 

0.066 

0.081 

0.40 

0.047 

0.067 

0.082 

0.45 

0.048 

0.068 

0.083 

0.50 

0.048 

0.070 

0.085 

0.900 

0.950 

0.975 

0.990 

0.084 

0.102 

0.119 

0.143 

0.085 

0.103 

0.121 

0.144 

0.086 

0.104 

0.122 

0.146 

0.087 

0.105 

0.123 

0.148 

0.088 

0.107 

0.125 

0.150 

0.089 

0.108 

0.127 

0.152 

0.091 

0. 110 

0.129 

0.154 

0.092 

0.111 

0.131 

0.157 

0.094 

0.113 

0.133 

0.159 

0.095 

0.115 

0.136 

0.162 

0.097 

0.118 

0.138 

0.166 

3* 


0.995 


0.187 

0.188 

0.189 

0.192 

0.196 

0.201 

0.209 

0.220 

0.235 

0.253 

0.253 


0.995 


0.160 

0.162 

0.164 

0.166 

0.168 

0.171 

0.173 

0.176 

0.180 

0.183 

0.187 
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